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ABSTRACT 

Context. Evidence has grown over the past few years that the neutral phase of blue compact dwarf (BCD) galaxies may be metal- 
deficient as compared to the ionized gas of their H n regions. These results have strong implications for our understanding of the 
chemical evolution of galaxies, and it is essential to strengthen the method, as well as to find possible explanations. 
Aims. We present the analysis of the interstellar spectrum of Pox 36 with the Far Ultraviolet Spectroscopic Explorer (FUSE). Pox 36 
was selected because of the relatively low foreground gas content that makes it possible to detect absorption-lines weak enough that 
unseen components should not be saturated. 

Methods. Interstellar lines of H i, N i, O i, Si n, P n, Ar i, and Fe n are detected. Column densities are derived directly from the 
observed line profiles except for H i, whose lines are contaminated by stellar absorption, thus needing the stellar continuum to be 
removed. We used the TLUSTY models to remove the stellar continuum and isolate the interstellar component. The best fit indicates 
that the dominant stellar population is B0. The observed far-UV flux agrees with an equivalent number of ~ 300 B0 stars. The fit of 
the interstellar H i line gives a column density of io 20 3±0 4 cm" 2 . Chemical abundances were then computed from the column densities 
using the dominant ionization stage in the neutral gas. Our abundances are compared to those measured from emission-line spectra in 
the optical, probing the ionized gas of the H n regions. 

Results. Our results suggest that the neutral gas of Pox 36 is metal-deficient by a factor ~ 7 as compared to the ionized gas, and they 
agree with a metallicity of ss 1/35 Z Q . Elemental depletion is not problematic because of the low dust content along the selected lines 
of sight. In contrast, the ionized gas shows a clear depletion pattern, with iron being strongly depleted. 

Conclusions. The abundance discontinuity between the neutral and ionized phases implies that most of the metals released by con- 
secutive star-formation episodes mixes with the H i gas. The volume extent of the enrichment is so large that the metallicity of the 
neutral gas increases only slightly. The star-forming regions could be enriched only by a small fraction (~ 1%), but it would greatly 
enhance its metallicity. Our results are compared to those of other BCDs. We confirm the overall underabundance of metals in their 
neutral gas, with perhaps only the lowest metallicity BCDs showing no discontinuity. 
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1. Introduction 

In the chemical downsizing scenario, in which massive galaxies 
are the first to form stars at a high rate (e.g., Cen & Ostriker 
1999), dwarf galaxies can remain chemically unevolved until 
later times. It is thus theoretically possible that some dwarf 
galaxies in the nearby Universe still contain pristine gas. The 
first observational constraint was provided by studies of chemi- 
cal abundances in the ionized gas of blue compact dwarf galax- 
ies (BCDs). These studies yield a fundamental result: a metal- 
licity floor appears to exist, with no BCD having an ionized 
gas oxygen abundance lower than 12 + log(0/H) ~ 7.0, or 
Z /5CO|. For more than three decades, I Zw 18 discovered by 
Zwicky (1966) and first studied by Sargent & Searle (1970) held 
the record as the most metal-deficient galaxy known, with an 
oxygen abundance 12 + log(0/H) = 7.17 ± 0.01 in its NW 
component and 7.22 ± 0.02 in its SE component (Thuan & 
Izotov 2005). Only very recently has IZw 18 been displaced 
by two BCDs that are more metal-deficient, SBS 0335-052W 
with 12 + log(Q/H) = 7.12 ± 0.03 (Izotov et al. 2005) and 



1 We use hereafter the solar oxygen abundance 12 + log(0/H) = 8.66 
of Asplund et al. (2005) 



DDO68 with 12 + log(0/H) = 7.14 + 0.03 (Izotov & Thuan 
2007). The difficulty of finding more metal-poor objects was 
originally explained by the local mixing of the amount of met- 
als produced and released in a single starburst episode (Kunth & 
Sargent 1986). This hypothesis has been since challenged by ob- 
servations of uniform abundances in disconnected star-forming 
regions within a single dwarf galaxy (e.g., Skillman & Kennicut 
1993; Kobulnicky & Skillman 1997; Noeske et al. 2000). 

Tenorio-Tagle (1996) proposed an attractive scenario in 
which supernovae (SNe) products are driven to high galactic lat- 
itudes in a hot phase (~ 10 6 K) before they cool down, con- 
dense as molecular droplets, and rain back on the disk. The 
complete cycle (before mixing) would last several ~ lOOMyr. 
This scenario naturally explains the uniform abundances found 
in dwarf galaxies, while ruling out local enrichment within the 
star-forming regions. It must be noted that models of Rieschick 
& Hensler (2001; 2003), based on Tenorio-Tagle's scenario, pre- 
dict that only ~ 75% of the SNe products go through the outer 
galactic cycle as hot gas, while the remaining 25% evolve in an 
inner local cycle and can enrich the star-forming regions in only 
lOMyr. A few percent could eventually leave the gravitational 
field of the galaxy. On the other hand, Chandra X-ray observa- 
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tions of SBS 0335-052W, SBS 0335-052E, and IZw 18 do not 
reveal any sign of a hot gas breaking out from the stellar body 
implying that the galactic fountain scenario may not be valid 
for BCDs (Thuan et al. 2004). In fact, Recchi et al. (2001; 2004) 
modeled star-formation in gas-rich dwarf galaxies and found that 
the majority of newly produced metals resides in a cold phase 
(< 2 x 10 4 K) after only a few tens of Myr, supporting the com- 
mon instantaneous mixing hypothesis. 

To understand the metal enrichment of galaxies, it is neces- 
sary to study all the gaseous phases involved in the gas mixing 
cycle. X-ray observations of starburst galaxies led to the con- 
clusion that the hot gas is probably enriched by SNe products 
(e.g., Martin et al. 2002). The following step in the mixing cycle 
is poorly constrained observationally. Do metals mix with the 
cold halo and thus enriching the H i region? Do they eventually 
condense as molecular droplets and settle back onto the stellar 
body? The situation might be different in BCDs if galactic foun- 
tains do not form. However, it seems that winds from massive 
stars and SNe are still able to carve the ISM to form superbub- 
bles in which metals are injected at high temperatures (Thuan et 
al. 2004). Mixing could then occur at the interfaces between the 
hot and cold gas. Thus the study of the metal content in the neu- 
tral gas of dwarf galaxies provides us with essential constraints 
for their chemical evolution models. 

The neutral phase usually consists of at least two distinct 
components (Field et al. 1969), with a warm medium (WNM; 
T w ~ 10 4 K, n w ~ few tenths crrT 3 ) showing broad hyperfine 
21cm line emission, and a cold medium (CNM, T c ~ 100K, 
n c ~ few tens cirT 3 ) showing narrow 21cm emission and ab- 
sorption. It must be noted that the distinction between CNM and 
WNM becomes indistinct in galaxies with low-pressure ISM, 
such as dwarfs (Lo et al. 1993). Neutral metals can be ob- 
served through resonant lines in the far-ultraviolet (FUV). BCDs 
are ideal targets because they display large amounts of H i gas 
(Thuan & Martin 1981), and because the massive stars provide 
a strong FUV continuum. Kunth et al. (1994) used the GHRS in- 
strument on board the Hubble Space Telescope to observe inter- 
stellar absorption lines along the lines of sight toward the many 
massive stars in IZwl8. The authors detected the O i line at 
1302.2 A and, together with the neutral hydrogen content mea- 
sured from the 21 cm line, estimated the metallicity of the neu- 
tral envelope to be at least 10 times lower than that of the ionized 
gas in the H n regions. Later, Thuan & Izotov (1997) obtained 
a GHRS spectrum for SBS 0335-052E and found the metallic- 
ity of the neutral gas to be similar to that of the ionized gas, 
about 2% solar. These results remained however inconclusive 
because of possible saturation effects of the O i absorption line 
(see Pettini & Lipman 1995). 

A new step forward was achieved with the launch of the Far 
Ultraviolet Spectroscopic Explorer (FUSE; Moos et al. 2000), 
which allows the observation of absorption-lines of H i together 
with many metallic species such as N i, O i, Si n, P n, Ar i, 
and Fe n. The spectral resolution of FUSE (R ~ 20 000) and 
its sensitivity made it possible to detect numerous absorption- 
lines, some of them apparently not saturated. The neutral gas 
chemical composition was derived in several BCDs, spanning a 
wide range in ionized gas metallicity, from 1/50 to 1/3 solar, with 
IZw 18 (Aloisi et al. 2003; Lecavelier des Etangs et al. 2004), 
SBS 0335-052 (Thuan et al. 2005), IZw 36 (Lebouteiller et al. 

2004) , Mark 59 (Thuan et al. 2002), NGC625 (Cannon et al. 

2005) , and NGC 1705 (Heckman et al. 2001). These investiga- 
tions showed that the neutral gas of BCDs is not pristine, but 
that it has already been enriched with metals up to an amount of 
> 1/50 Z . The second most important result is that the metal- 



licity of the neutral gas is generally lower than that of the ion- 
ized gas - except for the two lowest-metallicity BCDs, IZw 18 
and SBS 0335-052E - implying that the neutral phase has been 
probably less processed. 

The interpretation of these remarkable results has been lim- 
ited so far by the lack of knowledge about the origin of the neu- 
tral gas in BCDs whose absorption is detected in the FUV. Direct 
comparison with the H i gas seen via the 21 cm emission is ham- 
pered by a bias in selecting out relatively dust-free FUV lines 
of sight toward the massive stars. Even more importantly, FUV 
observations only probe gas in front of stellar clusters while ra- 
dio observations in emission probe the whole system. Many un- 
certainties also remain concerning the method of deriving col- 
umn densities from absorption-lines in complex systems. Most 
FUSE studies of BCDs assumed a single homogeneous line of 
sight with the exception of Mark 59 (Thuan et al. 2002; see also 
Lebouteiller et al. 2006 for the giant H n region NGC 604). A 
more realistic approach would be to consider multiple lines of 
sight toward massive stars, each line of sight intersecting clouds 
with possibly different physical conditions (such as turbulent ve- 
locity b, radial velocity v) and chemical properties. If individual 
absorption components are unresolved, it is almost impossible to 
tell whether they are saturated, in which case the column density 
inferred from the global absorption line can be severely underes- 
timated (see e.g., Lebouteiller et al. 2006). Following the sugges- 
tion of Kunth & Sargent (1986), Bowen et al. (2005) got around 
the problem by deriving the abundances in the neutral gas of the 
dwarf spiral galaxy SBS 1543+593 along a quasar line of sight. 
The authors found that the abundances agree with those in the 
ionized gas of the brightest H n region in the arms, contrasting 
with the FUSE results in BCDs. As far as chemical models are 
concerned, Recchi et al. (2004) examined the case of IZw 18 for 
which they found that abundances in the H i medium (defined by 
a temperature lower than 7 000 K) are similar to those in the H n 
gas. Hence, although some physics is still missing in the models 
(radiative transfer, clumpy ISM, dust phase), it seems that FUSE 
results still cannot be reproduced. 

In this paper, we look anew at the problem of how the 
chemical compositions of the ionized and the neutral gas in 
BCDs compare by selecting an object with an extremely low 
H i column density. This object has some absorption lines that 
are weak enough to be safely considered as unsaturated. The 
object we analyze is Pox 36 (also known as ESO572-G034, 
I SZ63, HIPASSJ1 158-19b, or IRAS 1 1564-1844), an emission- 
line galaxy first studied by Rodgers et al. (1978) and then re- 
discovered by Kunth et al. (1981). It is characterized by strong 
optical emission-lines superimposed on a very blue continuum. 
Its spectral characteristics and morphological properties fit the 
definition of a BCD galaxy (see e.g., Thuan & Martin 1981, 
Kunth & Ostlin 2000). Its known observational properties are 
listed in Tablefl] In particular, chemical abundances in the ion- 
ized gas of the H n regions within Pox 36 have been derived 
by several independent studies (e.g., Kunth & Sargent 1983; 
Izotov & Thuan 2004). The ionized gas has an oxygen abun- 
dance 12 + log (O/H) m 8.05 which corresponds to Z Q /5. With a 
H i velocity of 1 1 14 km s _I (Meyer et al. 2004), a Hubble con- 
stant of 73 kms _I Mpc~', its cosmology-corrected distance is 
20Mpc. 

The FUSE observations are presented in Sect.|2] and an 
overview of the FUSE spectrum is given in Sect. [3] We infer the 
column densities of H i and of the metals from the line profile 
in Sect. [4] Chemical abundances are derived and discussed in 
Sect. [5] Finally, we study the properties of the neutral gas in the 
current FUSE BCD sample to discuss general trends in Sect.|7] 
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2. Observations 

Pox 36 was observed by the FUSE telescope on May 11 2001 
(observation B0220201). All eight observing channels were used 
with the LWRS aperture (30" x 30"). At the distance of Pox 36 
(see Table[TJ, 1" corresponds to m 100 pc. The total exposure 
time is 25 470 s. Data were taken in "Time Tag" observing mode. 
The brightest blue regions are included in the aperture as well as 
the fainter extended emission seen toward the north-east (Fig.Q]). 
In the blue band, the combined emission from the bright knots is 
on the same order as the emission from the rest of the body. We 
expect shorter wavelength light such as the FUV emission to be 
dominated by the knots where most of the massive stars lie. 

Data were processed by the CaIFUSE0 pipeline (version 
3.2.1). We aligned the individual exposures on the wavelength 
grid and co-added them to obtain one spectrum per observing 
channel. The final spectrum is presented in Fig. [2] in which ob- 
servations from different channels have been co-added for dis- 
play purposes. The signal-to-noise ratio per resolution element 
(20 km s -1 ) is *T2. 

The intrinsic spectroscopic resolving power of FUSE is close 
to R — 20 000 corresponding to a line spread function (LSF) 
width of ~ 15kms _1 (e.g., Sahnow et al. 2000; Hebrard et al. 
2002; Rruk et al. 2002). The effective resolution depends on the 
quality of the individual exposure alignment. In particular, the 
alignment is difficult to assess for spectra in the SiC channels be- 
cause of the smaller effective area as compared to the LiF chan- 
nels, resulting in a lower signal-to-noise ratio. Misalignments 
reduce significantly the spectral resolution of the final spectrum 
so that we estimate the LSF width to be as large as ~ 25 km s _1 . 

The spatial extent of the source within the LWRS aper- 
ture can also result in a smearing of the spectrum (see e.g., 
Lebouteiller et al. 2004; 2006). Unfortunately, Pox 36 has not 
been imaged in the UV. The optical blue image of Pox 36 shows 
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Fig. 1. Blue band image of Pox 36 (Astro- Wise, ESO-LV 
Catalog). The separation between the 2 bright knots is given by 
the segment. The inclined square indicates the 30" x 30" LWRS 
aperture. 



two bright knots separated by m 6" (Fig.[TJ, corresponding to 
a projected distance of m 4" in the dispersion direction. Given 
the LWRS aperture size of 30" = 120kms _I , the separation 
between the 2 knots should result in a wavelength smearing of 
~ 16kms _1 at most. 

Since the line broadening also depends on several observa- 
tional parameters (such as the pointing jitter), it is difficult to de- 
rive a priori a precise estimate. We empirically determined the 
LSF width by using the separation between closed or blended 
absorption lines. The most likely value is ~ 30kms _1 for the 
LiF channels and ~ 37.5 km s~' for the SiC channels. This is 
roughly consistent with the quadratic sum of the intrinsic PSF 
and the smearing due to the spatial extent. This is the LSF width 
we use for line fitting (Sect. 14. It . In practice, our results depend 
little on the LSF width determination (see Sect. 14.1. Ti l. 



3. Overview of the FUSE spectrum 

3. 1 . Absorption systems 

The spectral continuum is provided by the UV-bright massive 
stars in the galaxy. Absorption lines from species along the line 
of sight are superimposed on the continuum. The line of sight 
intersects the ISM from the Milky Way and from Pox 36 itself. 
As noted by Wakker et al. (2003), the line of sight passes at only 
18' (a 100 kpc projected) from NGC4027, but there is no ab- 
sorption visible at the velocity of the latter (» 1460 km s~'). The 
redshift of Pox 36 makes it possible to separate easily the local 
absorption system from the one in the BCD. We identify ab- 
sorption lines from the Milky Way at nearly zero radial velocity. 
We barely detect Galactic H2, which is mostly due to the high 
Galactic latitude of Pox 36 (Table 1). 

The absorption system corresponding to the neutral ISM in 
Pox 36 is detected at a velocity of v„ ~ 1052 km s~' (1058 ± 
10kms _I from the metal lines, 1047 + 10km s" 1 from the H 1 
lines). This corresponds to a wavelength shift of ~ 3.5 - 4.0 A 
from the rest-frame. Only the strongest atomic lines are detected. 
No H2 line at the redshift of Pox 36 is seen. The velocity v„ is 
significantly smaller than the one derived from the 2 1 cm H 1 
line equal to 1113.9 ± 6.5 km s -1 (Meyer et al. 2004; Zwaan et 
al. 2004). However, it must be stressed that the regions probed 
in the UV and the radio are different in terms of extent (the radio 
beam size is 15.5' as compared to the 30" x 30" aperture of the 
LWRS) and of depth (because of dust extinction). The aperture 
(and depth) effect explanation is supported by the fact that v„ 
is in agreement with the value of 1065 ±20 km s _I derived by 
Izotov & Thuan (2004) from the optical emission lines of the 
ionized gas (Table[TJ). Taking at face the velocity of the neutral 
and ionized gas tracers, we conclude that there is no significant 
evidence that the neutral gas is outflowing from the star-forming 
regions. 

We find that the turbulent velocity of the neutral gas in 
Pox 36 is b„ = 33 ± 10 km s . Its value is mostly constrained by 
lines near saturation such as the O 1 988.77 A line (Table|2|. The 
quoted uncertainty on the b-value does not include possible sys- 
tematic errors on the LSF width. The corresponding FWHM is 
is 2b„ Vln2 w 55 ± 17kms _I , consistent with the FWHM mea- 
sured in radio (67.4 km s~', Table[TJ). However as said before, 
the region probed by the large radio beam (15.5') is not directly 
comparable to the gas probed in the FUV (corresponding to an 
extent of less than ~ 10", Fig.[TJ. It must be also stressed that 
our determination of b„ may not be physical if the absorption 
lines arise from a collection of interstellar clouds (see Sects. l3~2l 
and |4T2l . 
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Table 1. Observational properties of Pox 36. 



Property 


Value 


Comment 


a (J2000) 


Ilh58m57.3s 


HIPASS (Meyer et al. 2004) 


S (J2000) 


-19°01'41" 


HIPASS (Meyer et al. 2004) 


Galactic coordinates 


l in/ / a rt 1/11 

/ = 286.114, b = 42.123 


T TTD A C C /A X , 1 TAfl A \ 

HIPASS (Meyer et al. 2004) 


Foreground Galactic extinction (mag) 


t{B — V) = U.Uio 


C~U1~~~1 ~f ^1 /100Q\ 

Schlegel et al. (199o) 


Magnitude (mag) 


Bj - 14.27, / = 14.38 


HIPASS (Doyle et al. 2005) 




K s = 12.63 


Vanzi & Sauvage (2006) 


r /tt/j\ / — 1 \ 

L(H/?)(ergss L ) 


1 o i , , 1 r\ i 1 -' 

1.81 x 10 


From F(H/s) (lzotov & Thuan 2004) and D = 20 Mpc 


FWHM (H i) (kms -1 ) 


67.4 


HIPASS (Meyer et al. 2004) 


Peak H i flux (Jy) 


0.1197 


HIPASS (Meyer et al. 2004) 


T , . Ill 13 /T 1 — 1 \ 

Integrated H i flux (Jy km s ) 


8.6 


T 1 1 n A fin /"» n .1 '-y f\f\ a \ 

HIPASS (Meyer et al. 2004) 


M(H i) (Mo) 


8.1 x 10 8 


Derived from the H i flux with the distance=20 Mpc 


Radial velocity ( km s~' ) 


1114 ±2 


21 cm line, Theureau et al. (1998) 




1113.9 ± 6.5 


HIPASS (Meyer et al. 2004; Zwaan et al. 2004) 




1065 ± 20 


Optical lines (lzotov & Thuan 2004) 


Distance (Mpc) 


« 12 


Vigroux et al. (1987) with H = 73kms~ 1 




* 16 


Kunth& Joubert(1985) 




« 20 


Using v = llWkms- 1 , H = 73kms-», n„, = 0.27 


logO/H 


-3.95 ± 0.02 


lzotov & Thuan (2004) 




-3.86 ±0.03 


Kobulnicky & Skillman (1996) 




-3.81 ±0.04 


Campbell (1988; 1992) 




-3.82 ± 0.04 


Kunth & Sargent (1983) 


logN/H 


-5.53 ±0.03 


lzotov & Thuan (2004) 




-5.39 ±0.10 


Kobulnicky & Skillman (1996) 




-5.38 ±0.05 


Campbell (1992) 




-5.37 ± 0.05 


Kunth & Sargent (1983) 


log Ar/H 


-6.24 ± 0.03 


lzotov & Thuan (2004) 


log Fe/H 


-6.17 ±0.10 


lzotov & Thuan (2004) 



le-13 




1050 



1120 



1130 



1140 1150 
Wavelength (A) 



1160 



1170 



1180 



Fig. 2. Complete FUSE spectrum of Pox 36. The flux is plotted against the observed wavelength (3.5-4 A shift from the rest wave- 
length, Sect. 13. p . Data from overlapping channels were co-added for display purposes. The strong emission-lines at ~ 972 A and 
~ 1026 A correspond to terrestrial H i airglows. 
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5 



Specie 


<U (A) 


<U (A) 




N sst (cm- 2 t 


Comment 


Hi 


1025.722 


1029.389 


0.0527 




Lyman [} 


Hi 


937.803 


941.080 


0.0026 




Lyman e 


Ni 


953.415 


956.797 


0.0132 


5 x 10 14 




Ni 


953.655 


957.031 


0.0250 


3 x 10 14 




N i 


953.970 


957.343 


0.0350 


2x 10 14 




N i 


963.990 


967.040 


0.0148 


4 x 10 14 


Blended with P n ,1963.8 


N i 


964.626 


968.029 


0.0094 


7 x 10 14 


Barely detected 


N i 


965.041 


968.458 


0.0040 


2 x 10 15 


Not detected 


N i 


1134.165 


1138.069 


0.0134 


4 x 10 14 


Barely detected 


N i 


1134.415 


1 138.303 


0.0268 


2 x 10 14 




Ni 


1134.980 


1138.888 


0.0402 


8 x 10 13 




Oi 


988.773 


992.209 


0.0465 


1 x 10 14 


Saturated 


Oi 


1039.230 


1042.961 


0.0092 


5x 10 14 




Sin 


1020.699 


1024.358 


0.0164 


2x 10 14 




Pn 


963.800 


967.249 


1.4600 


2x 10'" 


Blended with N i ,1963.99 


Pn 


1152.818 


1156.906 


0.2450 


2 x 10 13 


Not detected 


An 


1048.220 


1051.957 


0.2630 


9x 10 1 - 




Fe n 


1144.938 


1148.976 


0.1060 


3 x 10" 





" Atomic data are from Morton et al. (2003), except for Fe n Al 144.938 for which we take Howk et al. (2000). 
b N sat is the maximum column density before saturation of a central component with b = 2kms~' (Sec. l4.1.2l . 



3.2. Neutral ISM components 

The lines of sight observed by FUSE are optically thick to the 
Lyman continuum of the stellar clusters but the UV radiation is 
still able to reach us. This implies that the amount of foreground 
dust is either very small and/or it is clumpy, as expected for 
metal-poor objects (e.g., Lisenfeld & Ferrara 1998; Hirashita et 
al. 2002). Although FUSE observations are sensitive to both the 
CNM and the WNM as defined in McKee & Ostriker (1977) and 
Wolfire et al. (1995), the situation might be complex in BCDs 
since the relatively low ISM pressure could lead to a single H i 
phase, except in regions affected by stellar winds and supernovas 
shocks. Unfortunately, the spectral resolution of FUSE (Sect.|2]l 
is too small to isolate warm and cold neutral phases. 

In principle, the presence of molecular absorption could al- 
low us to constrain the physical conditions along the observed 
lines of sight. No diffuse H2 has been detected so far along lines 
of sight toward massive stars in BCDs (e.g., Vidal-Madjar et al. 
2000). The reasons are intrinsic to the H 1 phase since H2 has 
been detected in emission in the warm ISM of these objects (e.g., 
Vanzi et al. 2005; Thuan et al. 2005). Warm H2 has been since 
observed in several other BCDs by the Infrared Spectrograph 
onboard Spitzer (Lebouteiller et al. in preparation). The absence 
of H2 in the diffuse ISM of BCDs is due to the combined ef- 
fects of a low H 1 density in the neutral gas envelope, a large UV 
flux that destroys H2 molecules and a low metallicity that makes 
grains on which to form H2 molecules scarce (Vidal-Madjar et 
al. 2000). Diffuse H2 is detected in relatively more metal-rich 
star-forming objects, but the molecular fraction is very low, on 
the order of 10~ 4 (Hoopes et al. 2004). A small velocity offset 
(~ 10 km s~') between H2 lines and atomic lines has been found 
in NGC 625 (Cannon et al. 2005) and toward the giant H n re- 
gion NGC604 (Lebouteiller et al. 2006). Although the offset is 
small, this might indicate that atomic absorption-lines preferen- 
tially arise in regions where H2 is not able to form. Hence, it 
seems reasonable to assume that the neutral phase probed by 
FUV absorption-lines is rather warm and diffuse since dense 
clouds should show H2 absorption. In the following, we thus 
assume that the interstellar absorption-lines arise from a single 



phase dominated by turbulent motions. Since the filling factor is 
unknown, this phase could be a uniform H 1 envelope or a col- 
lection of H 1 clouds. 

3.3. Stellar contamination 

The FUV spectral continuum of most of the BCDs and giant H 11 
regions studied so far with FUSE is dominated by O stars. This 
is indicative of a starburst age younger than ~ lOMyr (Robert 
et al. 2003). O stars have relatively "clean" spectra, the most 
prominent stellar absorption lines being due to He 11, O 11- m, 
and N 11- m in the optical, and to O vi, N v, P v, S iv, C iv and 
Si iv in the FUV. Because of the high temperature, hydrogen 
is mostly ionized, and only weak H 1 photospheric lines can be 
observed which contribute to the - already saturated - core of 
the interstellar H 1 absorption line. Hence FUV spectra toward 
young starbursts allow a precise determination of the interstel- 
lar H 1 content. The only exception concerns extremely young 
starburst episodes showing a prominent O vi P-Cygni profile, 
which complicates a reliable profile fitting of the wings of Ly/3 
(Lebouteiller et al. 2006). 

Older starbursts (> lOMyr) are dominated by B stars. Such 
stars have particularly strong He 1 lines in the optical, and over- 
all a more complex FUV spectrum because the ionization degree 
is lower than in O stars. For this reason, photospheric H 1 lines 
becomes prominent (Valls-Gabaud 1993). As an illustration, the 
H 1 lines from the Lyman series in the BCD Mark 59 show "V- 
shaped" wings typical from B stars photospheres (Thuan et al. 
2002). Pox 36 shows a similar pattern, but while Thuan et al. 
(2002) got around the problem by considering an artificial con- 
tinuum to fit the profile of the H 1 lines in Mark 59, we decided 
to model the stellar absorption in Pox 36 (Sect. |4.2| i. 
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Table 3. Column densities. Errors are given at 2 <x. 



Specie 



(cm" 2 ) 



N b 

1 v norn^ 

(cm" 2 ) 



N i 
Oi 
Sin 
Pn 
Ar i 
Fe n 



14 20 +012 
i c co+0.12 

15.06^1 
12.78«| 
13.18tg 



13 96+ 022 
15 52 + < i "> 

i5.ooi;5 

12.94 +0 16 
13 10 + ^' 
13 97+° 13 



" N obs is the column density calculated using local spectral continuum 
from the observation. 

* Mmrm is the column density calculated using the spectrum corrected 
from stellar absorption (Sect. l4.2] l. 



4. Column densities 

In order to infer chemical abundances, we derived column den- 
sities based on the line profiles. We adjusted the profiles of the 
metallic lines using the observed spectrum (Sect. 14. It . The H i 
column density was calculated from the Ly/3 line in the spectrum 
corrected for stellar absorption (Sect. 14.2] ). 



4.1. Metals 

4.1 .1 . Single component analysis 

For the profile fitting, we first made the usual assumption of a 
single line of sight intersecting a medium with uniform proper- 
ties (namely with constant radial and turbulent velocities, and 
column densities). The actual situation can be more complex, 
with multiple lines of sight toward the many massive stars and 
H n regions within the galaxy (Sect. l4.l31 l. 

The method to derive column densities from FUSE line pro- 
files is discussed in Hebrard et al. (2002). Profiles are adjusted 
using the Owens routine (Lemoine et al. 2002) through a mini- 
mization of the x 2 between the model and the observations. All 
the lines could be fitted with Voigt profiles. We considered that 
all the neutral species within Pox 36 share the same turbulent 
and radial velocities. Those parameters as well as the column 
density of each species are free to vary. The local continuum 
around the lines of interest was approximated by a polynomial 
of degree to 3. The continuum level and shape are also free 
parameters. The LSF width (Sect.[2jl has little influence on the 
column density determination. This is principally due to the fact 
that the lines providing constraints on the column density are 
weak. The column density determination of all the species con- 
sidered in this study varies by less than 0. 1 dex when forcing the 
LSF width to vary between 20 km s and 45 km s _1 . Other LSF 
widths result in unsatisfactory fits based on the x 2 value. Lines 
were fitted simultaneously for the different observing channels 
to benefit from the data redundancy allowed by FUSE. The fit- 
ted lines are those in Table|2] Final fits are presented in Figs. [3] 
|4j and[5]while derived column densities are given in Table[3] For 
consistency, we also derived column densities of metals using 
the spectrum corrected for stellar absorption (Sect. 14. 21 



4.1.2. Presence of saturated unresolved components 

In addition to the simple case of a single line of sight intersect- 
ing a neutral cloud with uniform properties (Sect. |4. lTTb , we also 



consider the case of a collection of diffuse clouds producing mul- 
tiple unresolved absorption components. Such clouds can have 
different physical conditions and chemical properties. It is com- 
monly accepted that a complex absorption structure along a sin- 
gle line of sight usually does not introduce systematic errors on 
the column density determination (Jenkins 1996). However, the 
presence of multiple lines of sight inevitably complicates the 
analysis. The most obvious effect is the wavelength smearing 
due to the distribution of the lines of sight across the aperture 
(Sect.lUl. 

The unusually weak interstellar lines in Pox 36 may indicate 
that unresolved absorption components - if present - are not 
saturated. In such a case, the equivalent width of the integrated 
absorption line is simply the sum of the equivalent widths of 
the individual components, although the contribution from lines 
of sight toward stars with different UV brightnesses might not 
scale linearly (Labasque & Boisse in preparation). We check for 
the presence of possible saturated components by following the 
method described in Lebouteiller et al. (2006). We derive the 
maximum column density the narrow component (b — 2 km s _1 ) 
in the center of a line can have before becoming saturated. The 
resulting values are given in Tableland should be compared to 
the column densities inferred from the global profiles (Table[3j. 
It shows that O i and Si n lines could be strongly saturated if 
narrow components are present, that the Fe n line could be satu- 
rated, but that N i, P n, and Ar i lines are only close to saturation. 
Note that this represents only the worst-case scenario since in- 
tervening clouds might have in fact larger velocity dispersion. 
Nevertheless, these results emphasize the reliability of the N i, 
P n, and Ar i column density determinations. 



4.2. Hydrogen 

The interstellar H i lines are strongly contaminated by absorption 
from stellar atmospheres (Sect. l3.3l l. This can be best seen for 
the high-order lines of the Lyman series in Fig.|2]which display 
typical "V-wings". The derivation of the interstellar H i column 
densities is done in three steps: 



estimate the stellar H i absorption using the high-order 
Lyman lines (e.g., Lye) whose wings are purely stellar in ori- 
gin, 

correct Ly/3 (whose wings are partly stellar and interstellar) 

by subtracting the stellar H i absorption, 

estimate the interstellar H i profile using the corrected Ly/3 

line. 



We used the TLUSTY FUV absorption model spectra for O 
and B stars of Lanz & Hubeny (2003, 2007) to reproduce the 
spectral continuum of Pox 36. We chose to model the stellar 
spectrum of Pox 36 with a single stellar population, which we 
will show in Sect. l4.2.Tl is a safe assumption. TLUSTY mod- 
els are available for several temperatures and metallicities. We 
considered only main-sequence stars (log g = 4.0) and have in- 
cluded the effects of stellar rotation. The average rotational ve- 
locities have been taken from Kippenhahn & Weigert (1990) and 
Lang (1992). Finally, it must be noted that the velocity disper- 
sion of stars within Pox 36, or between individual H n regions, is 
negligible before the stellar rotation velocity. To compare with 
the observed spectra, we normalized the latter using a spline fit 
for the continuum, and the synthetic spectra using the featureless 
continuum provided by TLUSTY. 
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Fig. 6. The Lye line profile gives a strong constraint on the pho- 
tosphere temperature of the dominant stellar population. Models 
are drawn for a metallicity of 1/1 000 Z Q , but the profile of stellar 
H i lines is unchanged for varying metallicities. 



4.2.1. Stellar temperature 

The wings of the stellar H i lines are strongly dependent on tem- 
perature. The best constraint is given by the high-order lines of 
the Lyman series since the interstellar contribution in these lines 
is saturated and they have no damping wings. Moreover, H i lines 
are not affected by metallicity so that they are not prone to the 
degeneracy between temperature and metallicity. 

Given the redshift of Pox 36 and the resulting contamina- 
tion by metallic lines from the Milky Way (see full spectrum in 
Fig.[2]i, Lye is the cleanest H i line for use to constrain the stel- 
lar temperature. Its line profile gives a temperature of 30 000 K 
(Fig.[6]i, corresponding to a B0 stellar population. The uncer- 
tainty on this temperature determination is smaller than the steps 
allowed by the models (< 2 500K). Are other stellar types 
present in large quantities in Pox 36 ? Hotter stars would con- 
tribute narrow H i absorption lines, not modifying appreciably 
the "V-shaped" wings of the high-order Lyman lines. As for the 
wings of the low-order Ly/3 line, they are dominated by inter- 
stellar absorption. Cooler stars would produce more prominent 
"V-wings" which are not seen in the profiles of the high-order 
Lyman lines. We conclude that the stellar populations responsi- 
ble for the stellar H i profile are mainly composed of B0 stars, 
with a possible small amount of hotter stars. 

Based on the surface flux given by the models and assuming 
a distance of 20 Mpc (Table[T|i, the observed FUV flux of Pox 36 
implies that the equivalent number of B0 stars is ~ 300. This is 
the number of unextincted foreground stars. 

4.2.2. Stellar metallicity 

Lines from ionized metals (C m- iv, N in- v, Si iv- v, ...) arise 
in the stellar atmospheres of O and B stars. These lines can 
also contaminate the interstellar spectrum although their effect 
is much weaker as compared to the stellar H i lines. One expects 
that the metallicity of young stars to be similar to that of the H n 
gas (see e.g., Venn et al. 2001; Lee, Skillman & Venn 2006), 
being equal to ~ 1/5 Z (Tabled in the case of Pox 36. 

We computed TLUSTY models for a range of metallicities 
so they can be used to compare with observations. We use the 
prominent C in lines at 977.02 A and at 1175.6 A to derive an 
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Fig. 7. The profile of the C m lines at 977.02 A (top) and at 
1175.6 A (bottom) is strongly affected by the metallicity. The 
models are from bottom to top forZ = 1/2, 1/10, 1/50, 1/100, 
and 1/ 1000 Z . 

estimate of the stellar metallicity. The 977.02 A line is a res- 
onance line which can be either interstellar or stellar. It gives 
an upper limit to the stellar absorption. On the other hand, the 
1 175.6 A multiplet is not resonant and is purely stellar. Based on 
the two C in lines, we find that the observations agree best with 
a metallicity of 1/1000 Z (Fig. [7]). This model fails however 
to reproduce the stellar absorption lines between 11 10 A and 
1120 A (Fig. [8]l. Inversely, the 1/5 Z Q model reproduces some 
stellar lines while overestimating other lines (notably around 
~ 1145 A). 

Thus it seems difficult to estimate the stellar metallicity by 
comparing the observations with the stellar lines predicted by 
models. This is of little importance for the interstellar H i col- 
umn density determination, but it may prevent a reliable esti- 
mate of the continuum for measuring metal column densities. In 
practice, the stellar lines are broad enough that they don't affect 
significantly the profiles of the narrow interstellar lines. We de- 
cided to use the 1/5 Z Q TLUSTY model and to rederive column 
densities of the metals using the normalized stellar absorption- 
corrected spectrum. The results are given in the last column of 
Table[3] They are consistent within the errors with those derived 
from spectra without correction for stellar metallic absorption, 
with perhaps the exception of the N i column density. Part of the 
small difference is probably due to systematic errors introduced 
in the normalization of the observations. Because the two deter- 
minations agree relatively well, we decide to use from now on 
the metal column densities derived from the observations with- 
out correction for stellar metallic absorption (Sect. [5]). 

4.2.3. Interstellar H i column density 

The interstellar H i content is mostly constrained by the damp- 
ing wings of the Ly/3 line, since higher-order lines are satu- 
rated. In order to correct for the stellar absorption, we use the 
TLUSTY stellar model. We chose a metallicity of l/1000Z o 
instead of 1 /5 Z to avoid uncertainties on the metallic stellar 
lines which we already saw cannot reproduce well the observa- 
tions (Sect. 14.2^21 ). This is a safe choice for deriving interstellar 
H i column density since the H i line profiles are independent of 
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Fig. 8. The stellar models for Z = 1/1000Z G (solid line) and 1/5 Z (dashed line) are superimposed on the normalized observed 
spectrum (gray histogram). The most prominent stellar lines are labelled. 



metallicity. In practice, our results on the H i column density do 
not depend on the metallicity chosen for the stellar models. 

We thus built the interstellar spectrum by removing the stel- 
lar contribution. In Fig. [9] the observations are compared to syn- 
thetic profiles taking into account stellar and interstellar absorp- 
tions. More general cases with various stellar temperatures and 
lines other than Ly/3 can be found in the appendix [A] The best 
model agrees with log (H i) * 20.3. The final fit, after removal 
of the stellar absorption is shown in Fig.[l0] The most likely col- 
umn density is log N (Hi) = 20.28 ±0.06. The uncertainty is only 
statistical and does not include systematic errors on the stellar 
model. We studied the influence of the stellar model on the H i 
column density by considering models at 29 000 K and 32 500 K. 



The fits can be found in Fig.[l0] The H i column density for the 
29 000K models is 20.18 ± 0.08. For such a temperature, even 
Ly/Ts profile is almost independent of the interstellar contribu- 
tion. In fact, the interstellar H i profile for 28 000 K is already 
saturated without any damping wings, and the H i column den- 
sity is not constrained (log (H i) < 20.1). For 32 500 K, the H i 
column density is 20.66 ± 0.06. In that case, the higher temper- 
ature results in narrow stellar lines, and the interstellar contribu- 
tion, and damping wings are clearly seen. Since a temperature 
of 32 500 K is clearly too high to reproduce correctly high-order 
Lyman lines (Sect. |4.2TT >, we consider that the corresponding H i 
column density should be regarded as a conservative upper limit. 
It must be noted that the fit quality (given by the x 1 ) is almost 
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Fig. 9. The observed Ly/3 profile is compared to synthetic pro- 
files combining the stellar model (30 000 K) and interstellar H i 
profiles with various column densities. 



identical for the 3 temperatures we tested. The reason is that 
Ly/3 is not a strong constraint on the stellar temperature con- 
trary to Lye (Sect. l4.2TTT >. We use hereafter the interstellar H i 
column density logN (H i) = 20.3 ± 0.4 with conservative error 
bars which include uncertainties on the best stellar model, on the 
stellar model itself, and on the interstellar profile fitting. 

Is this determination in agreement with radio observations? 
Assuming the neutral gas to be optically thin and a distance of 
20 Mpc, the integrated H i radio flux (TableQ} gives a mass of 
M(HI) = 8.1 x 10 8 M Q . If the H i distribution is uniform, the col- 
umn density we derived implies that the H i diameter is 13kpc. 
Based on the image of Fig.Q] the angular diameter of the galaxy 
in the optical is close to 30", which at the distance of Pox 36 
corresponds to » 3 kpc, i.e., a factor * 4 smaller than the H i 
extent. This is consistent with H i observations of BCDs (e.g., 
Thuan & Martin 1981). We conclude that the H i column density 
measured in the FUV is in agreement with radio observations. 



5. Chemical abundances 

In this section we derive the chemical abundances from the col- 
umn densities. Abundance determinations are mostly affected by 
ionization corrections and by element depletion on dust grains. 

5.1. Ionization structure 

The abundance of an element is usually estimated from its pri- 
mary ionization state and by applying ionization correction fac- 
tors. We expect that all elements with ionization potentials larger 
than that of hydrogen (13.6eV) are in a neutral state in the H i 
gas. This is the case for N, O, and Ar. On the other hand, Si, P, 
and Fe are mostly found as single-charged ions with negligible 
fractions of neutral atoms. 

Some nitrogen and argon could be ionized in places where 
hydrogen is neutral. This is a consequence of the large photo- 
ionization cross-section of N i and Ar i together with the pene- 
tration of UV photons in low-density partly-ionized regions of 
the ISM (Sofia & Jenkins 1998). We expect the ionization cor- 
rection factor for N i to be much smaller than for Ar r. Our results 
suggest in fact that they are both negligible (Sect.[6]l. 
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Fig. 10. Profile fitting of LyyS in Pox 36 at the observed wave- 
length of ss 1029.4 A. The spectra are corrected for stellar ab- 
sorption using models at 29 000 K (top), 30 000 K (middle), and 
32 500 K (bottom). The corresponding line in the Milky Way can 
be seen at the nominal wavelength of ss 1025.72 A. 



Inversely, ionization corrections for Si, P, and Fe may be 
needed if the singly-ionized ions exist in the ionized gas of H n 
regions and contribute to the observed FUV absorption lines. 
For the giant H n region NGC604, Lebouteiller et al. (2006) 
found that only Si n and P n require significant ionization cor- 
rections (0.34 dex and 0.15 dex resp.; see also Aloisi et al. 2003 
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Table 4. Chemical abundances. [X/H] is defined as log (X/H) 
- log (X/H) , where log (X/H) is the solar abundance. Solar 
abundances are from Asplund et al. (2005). 



Element (tracer) log(X/H) [X/H]£ [X/H]£ 



N (N i) 


-6.10 


±0.42 


-1.88 ± 


0.42 


-1.31 


±0.06 


0(0 i) 


-4.72 


±0.41 


-1.38 ± 


0.41 


-0.61 


±0.05 


Si (Si n) 


-5.24 


±0.47 


-0.75 ± 


0.47 






P(Pn) 


-7.52 


±0.47 


-0.88 ± 


0.47 






Ar (Ar i) 


-7.12 


±0.47 


-1.30 ± 


0.47 


-0.42 


±0.06 


Fe (Fe n) 


-6.32 


±0.42 


-1.77 ± 


0.42 


-1.62 


±0.11 



" Abundance in the neutral gas (this study). 

b Abundance in the ionized gas (Izotov & Thuan 2004). 



in I Zw 18). Thuan et al. (2002) estimated the amount of Fe n in 
the ionized gas of Mark 59 to be 20% of the total iron abundance 
derived from the absorption-lines. Therefore, we could expect 
ionization corrections on the same order to exist in Pox 36. We 
neglect them in the following discussion because it has only little 
influence on our results. 

We finally use N i, O i, Si rx, P n, Ar i, and Fe n to derive 
chemical abundances (Table|4|. It must be reminded that the H i 
column density is larger than the metal column densities by sev- 
eral orders of magnitude, so that weak absorption components 
could contribute to the damped H i line Ly/? while the corre- 
sponding component in metal lines are below the detection limit 
(see Hebrard et al. 2003 for the effect in D/H determinations). 
This would lower the interstellar abundance determinations as 
compared to their true values (Sect.[6]l. 

5.2. Depletion on dust grains 

The refractory elements Si and Fe are likely to be locked on dust 
grains and thus depleted from the diffuse phase. To compare the 
abundances of the neutral and ionized gas, we need to take into 
account potential depletion effects. 

It is difficult to assess the depletion of elements in the ionized 
gas of Pox 36 since only one refractory element, Fe, has been 
observed in this phase. The use of iron as a probe of depletion 
is complicated by the fact that it might be released in the ISM 
on larger timescales. In the Milky Way, it is thought that about 
2/3 of the iron in the ISM was produced by SNe type la involv- 
ing low-mass stars. Rodriguez & Esteban (2004) and Izotov et 
al. (2006) observed that the [O/Fe] ratio increases with metal- 
licity in emission-line galaxies, being consistent with depletion 
of Fe on dust grains, and with iron production in SN type II in 
those galaxies. The [O/Fe] ratio in the ionized gas of Pox 36 is 
[O/Fe], = 1.12 + 0.12, i.e., much larger than the mean value of 
0.40 ± 0.14 in the ionized gas of BCDs (Izotov & Thuan 1999). 
Accordingly, this should imply a significant iron depletion on 
dust grains in this phase. 

The dust content in the neutral gas of BCDs is unfortu- 
nately largely unknown. FUV observations only select lines of 
sight with little dust content. Vidal-Madjar et al. (2000) have ex- 
plained the absence of diffuse H2 in the BCD IZwl8 by by the 
low abundance of dust grains (together with the low H 1 vol- 
ume density and the hard radiation field). To date, no diffuse H2 
has been detected in other BCDs, including Pox 36 (Sect. 13. It . 
Therefore, we expect Fe to be less depleted in the neutral gas 
of Pox 36 than in the ionized gas. This is the case: the [O/Fe],, 



in the neutral gas of Pox 36 is 0.39 + 0.16, much lower than the 
value in the ionized gas (1.12 + 0.12), but consistent with the 
mean value in the ionized gas of BCDs (Izotov & Thuan 1999). 
We conclude that iron is nearly not depleted on dust grains in the 
neutral gas of Pox 36. Furthermore, since silicon is usually less 
depleted than iron (e.g., Sofia, Cardelli & Savage 1994; Savage 
& Sembach 1996), our silicon abundance determination should 
also be little affected by dust depletion. 



6. Results 

Based on the metal abundances of Table|4] it is clear that the 
neutral gas in Pox 36 is not pristine: it has already been enriched 
with heavy elements. This is of course expected since Pox 36 is 
more metal-rich than IZw 18 and SBS 0335-052 where no pris- 
tine gas has been found either (Aloisi et al. 2003; Lecavelier et 
al. 2004; Thuan et al. 2005). The origin of metals in the neutral 
phase is however unclear. The presence of old stellar population 
in the great majority of BCDs (e.g., Aloisi et al. 1999; Ostlin 
2000) leads to the conclusion that previous enrichment could 
have occured in previous starburst episodes or through low-level 
star-formation over the galaxy's history. Comparison with the 
chemical composition of the ionized gas in the H 11 regions pro- 
vides an interesting diagnostic to understand the enrichment pro- 
cess of the neutral gas. 

6. 1 . Comparison with the ionized gas of the H n regions 

Table|4] shows that N is underabundant by a factor m 4 in the 
neutral gas as compared to the ionized gas, while both O and 
Ar are underabundant by a factor as 6 and « 8 respectively. The 
relatively good agreement between the deficiency factors of N, 
O, and Ar in the neutral gas indicates that the abundances of N 
and Ar are probably well determined (Sect. 15. It . The iron abun- 
dance appears to behave differently, being the same within the 
errors in the two gaseous phases. However, if we correct [O/Fe] 
in the ionized gas of Pox 36 for dust depletion (Sect. |5.2t , the 
iron abundance becomes: 

[Fe/H]* ~ [Fe/H],. + [O/Fe],. - 0.40, (1) 

giving [Fe/H] ; ~ -0.9, about 7 times higher than the value in 
the neutral gas. Thus Fe is underabundant in the neutral gas by 
the same factor as N, O, and Ar. 

Can this metal deficiency in the neutral gas be attributed to an 
overestimate of the H 1 column density? Considering the metal 
column densities of Table[3] the neutral gas would have the same 
metallicity as the ionized gas if the H 1 column density is on the 
order of logN (H 1) ss 19.6. This is inconsistent with the ob- 
servations since such a low column density does not reproduce 
the wings of Ly/3 (Fig.|£t. We conclude that the discontinuity be- 
tween the metallicity of the neutral and ionized gas is not driven 
by uncertainties on the H 1 column density determination. 

Oxygen and argon are the best metallicity tracers available in 
the ionized gas, and agree with a metallicity of « 1/5 Z Q in this 
phase. Considering these 2 elements, we find that the metallicity 
in the neutral gas is a; 1/35 Z , i.e., a factor ~ 7 below the value 
in the ionized gas. Our results for Pox 36 are similar to those 
obtained for other BCDs studied by FUSE (Sect.[7ji: although the 
neutral gas of Pox 36 is not pristine, it is less chemically evolved 
than the ionized gas in the H 11 regions. 



14 V. Lebouteiller et al.: Metal enrichment of the neutral gas of blue compact dwarf galaxies: the compelling case of Pox 36 



6.2. Pristine gas pockets 

It is possible that pockets of pure H i lie in the halo, while only 
a fraction of the H i has been enriched close to the starburst re- 
gion. This would dilute the chemical abundances inferred from 
integrated lines of sight by adding H i without any metallic coun- 
terpart. In Pox 36, the H i has a radial velocity slightly lower than 
that of the metals (Sect. 13. 11 1. The same result was observed in 
IZwl8 by Aloisi et al. (2003) where the H i has a velocity of 
753 ± 6kms~ I as compared to 765 ± lOkms -1 for the metals. 
If the offset is real (i.e., not driven by systematic uncertainties 
on the wavelength calibration), this suggests that part of the H i 
could lie indeed in regions without any metals. This would im- 
ply that the ISM of the galaxy is in fact not well mixed, with 
much lower abundances toward the external regions. The offset 
is however too small to assert unambiguously the presence of 
metal-free gas. 

6.3. Enrichment process 

We now calculate the amount of metals produced and released 
by the present starburst episode in order to understand the metal 
enrichment of the neutral gas. The stellar mass formed is de- 
rived from the UV luminosity, following Lequeux et al. (1981). 
We assume an initial mass function slope of -1.5 and a burst 
age of lOMyr (see Sect. l3.31 >. We find that the stellar mass is 
2.5 X 10 6 M Q . Yields were taken from Meynet & Maeder (2002), 
assuming a metallicity of 0.2 Z Q and a mass interval of [2 - 
120] M . We find that an oxygen mass of M(O) = 1 .25 x 10 4 M 
is eventually released by the starburst. Several scenarios are then 
possible, depending on the volume extent over which the disper- 
sal and the mixing occurs. 

We first assume that all the newly produced oxygen mixes 
with the ionized gas. Following Kunth & Sargent (1986), the 
mass of ionized gas is: 



M(HII) 



Nrrip 



(2) 



where N is the number of ionizing photons, m p is the proton 
mass, «, is the electronic and ionic density, and j3 is the recom- 
bination coefficient. From the Hfi luminosity (Table[T|), we find 
N - 4.4x 10 50 s _1 . We calculate /3 = 2.482 x 10~ 14 cm 3 s _1 from 
the Table of Hummer & Storey (1987) using the electron den- 
sity (60 cirT 3 ) and temperature (12 560 K) found by Izotov et al. 
(2004). The mass of ionized gas is then M(HII) = 2.5 x 10 5 M Q . 
If the SNe products mix only with the ionized gas, the metallic- 
ity of this phase would reach an additional oxygen enrichment 
of 12 + log(0/H) = 9.64. This unrealistic value suggests that 
heavy elements will mix with a larger mass of gas, at larger spa- 
tial scales. 

Next, we assume that the newly produced oxygen will mix 
uniformly with the H i gas. Using the H i mass from TableQ] we 
find that the neutral gas would be enriched by an additional 5% 
only (12 + log(0/H) = 5.98). Therefore, it can be seen that many 
starburst episodes (~ 20) are required for the neutral gas to reach 
the metallicity level we derived with FUSE. 

A possible alternative is that the mixing of new metals might 
not be uniform as the H i lying close to the stellar clusters could 
be more enriched than the outskirts which would remain metal- 
free. In order for a nearby H i gas shell to have the same metal- 
licity as that of the star-forming region (ionized gas), we find 
that the newly produced oxygen should mix with 7.0 x 10 7 M 
of H i, i.e., ~ 10% of the total H i in the galaxy. This can be 
seen as the upper limit of the H i volume that can be enriched 



to reasonable levels. We conclude that SNe products could mix 
with a small fraction of the H i gas in Pox 36 while not enriching 
most of it. This hypothesis implies that pockets of pristine gas 
could well be present in the galaxy. One way to prove this is to 
study absorption-lines toward several star-forming regions of a 
given galaxy, with significantly different foreground gas content. 
Another way is to analyze the scatter of the relation between 
abundances in the ionized gas and neutral gas for a sample of 
BCDs (SectEB. 

Toward a plausible scenario, it could thus be that the vast 
majority of the metals released during a starburst episode mixes 
with a large mass of H i, increasing only slightly its metallicity, 
while a small fraction mixes with the star-forming region itself 
(currently ionized), increasing greatly its metallicity. These met- 
als most likely come from previous episodes of star-formation, 
with the consequence that star-formation episodes have to be 
co-spatial. We consider the simple case in which the metallic- 
ity of the neutral phase increases by f„m/M(HI), where m is the 
mass of metals released by a starburst and f„ is the fraction of 
these metals to be mixed within the H i region. Similarly, the ion- 
ized gas increases its metallicity by fiin/M(HII). Over succes- 
sive starburst episodes, a metallicity offset between the ionized 
and neutral phases is achieved and maintained if the following 
condition is met: 



fn 



1 



1 + a 



M(HII) 
M(HI) 



(3) 



where a is a parameter depending on the metallicity offset. An 
offset of 1 dex - as it is observed in BCDs - implies a = 10. If, 
as in Pox 36, the H i mass is 1000 times larger than the ionized 
gas mass, the fraction of metals mixing with the H i phase should 
be * 99%, while only 1% mixes with the gas of the star-forming 
region. In conclusion, under the previous condition, the neutral 
gas can remain metal-poor, compared to the ionized gas, as long 
as the star formation keeps occurring at the same location. If an- 
other burst occurs at another location, one might expect to mea- 
sure low abundances, but the offset will reappear when the next 
starburst episode occurs and after mixing is achieved. 

7. Discussion 

7.1. Chemical enrichment of BCDs 

The FUSE BCD sample is now large enough to allow the study 
of general trends. Are the results concerning Pox 36 consistent 
with those of other BCDs? A total of 7 BCDs with various 
metallicities (1/50 to 1/3 Z ) have now been analyzed, in ad- 
dition to the high-metallicity giant H n region NGC 604 in the 
spiral galaxy M33. We plot in Fig.Q~T]the elemental abundance 
[X/HlJ in both the ionized and neutral gas of the BCDs. It can 
be seen that N, O, and Ar show identical trends, with systemat- 
ically lower abundances in the neutral ISM than in the ionized 
gas. Part of the dispersion of [O/H] is probably due to saturation 
effects which could not be avoided in some objects. Some of 
the dispersion of [N/H] and [Ar/H] is probably due to ionization 
corrections. Since the trend is similar for N, O, and Ar, the cor- 
rections are either of the same order for all objects (< 0.5 dex), 
or they are negligible, as for Pox 36 (Sect.[6]l. The abundance of 
iron in the neutral and ionized gas of BCDs must be interpreted 
with more care because of depletion on dust grains (Sect |5.2t . 
Figure[T_2] shows the relation between the Fe abundances in the 



3 [X/H] is defined as log (X/H) - log (X/H) , where log (X/H) G is the 
solar abundance from Asplund et al. (2005). 
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Fig. 11. The abundances of nitrogen (top), oxygen (middle), 
and argon (bottom) in the ionized gas and in the neutral gas 
are compared. The dashed line indicates the 1:1 ratio. Labels: 
(1*) IZwl8 (Lecavelier des Etangs et al. 2004), (1) IZwl8 



Fig. 12. Abundances of Fe in the neutral gas and in the ionized 
gas. Grey points correspond to the [Fe/H] in the ionized gas 
while black dots include a correction to account for depletion 
(Sect.iJl. See Fig.HUfor the labels. 



two gaseous phases, considering both the observed abundances 
and the ones corrected for depletion. It can be seen that corrected 
values agree more with the general trends observed for N, O, and 
Ar (Fig.rTTIi. It is seen in particular that iron depletion affects sig- 
nificantly the ionized gas abundances, especially in SBS 0335- 
052, Pox 36, and NGC 604. 

Examination of [Ar/H] and [N/H] in Fig.QT]reveals the pos- 
sible existence of a plateau for the lowest-metallicity BCDs 
(12+log(0/H), < 7.8 or < 1/8 Z ). The data on [O/H],, are 
also consistent with the existence of a metallicity plateau at 
12 + log(0/H)„ ~ 7.0 ([0/H]„ = -1.7), although the error bars 
are considerably larger. The presence of a plateau would imply 
a minimal enrichment of the neutral ISM. It could be that the 
H i reservoir has this minimal metallicity before star-formation 
ever occurred. Interestingly, Telfer et al. (2002) found that the 
intergalactic medium in the redshift range 1.6 < z < 2.9 has a 
metallicity 12 + log(0/H) between 6.7 and 7.6, which includes 
the floor metallicity observed in BCDs. Hence it appears that 
the neutral and ionized ISM of the lowest metallicity galaxies 
could well be characterized by a default enrichment rather than 
by previous star-formation episodes. In opposition, the higher- 
metallicity objects would be enriched by many starburst episode. 
Further observations and analysis of BCDs with metallicities 
lower than ~ 1 /8 Z will be needed to confirm the existence 
of this plateau. 

At higher metallicities > 1 /5 Z Q , the abundances of the neu- 
tral gas and ionized gas seem to be well correlated. This is the 
case for Mark59, Pox 36, NGC 625, and NGC 1705. The giant 
H ii region NGC 604 appears to extend the trend to nearly so- 
lar metallicities. The neutral gas abundances are systematically 
lower than those of the ionized gas. Nitrogen, oxygen, and argon 
in the neutral gas are underabundant by a factor ~ 0.9 dex, i.e., 
nearly 10, as compared to the ionized gas. 

If we dismiss the idea of a plateau and consider instead that 
the points in Fig.[TT]show a positive correlation with significant 
scatter from the lowest to the highest metallicity objects, then 
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the situation could be consistent with a dilution by metal-free 
gas along the lines of sight. The scatter would then be due to 
the dilution factor being different for each object. However one 
would then expect objects with the largest H i column density to 
display the largest underabundance in the neutral gas. We stud- 
ied the dependences of the underabundance factor with the H i 
column density but found no significant correlation. 

7.2. Interpretation 

The lower metallicity observed in the H i regions of BCDs as 
compared to the ionized gas of their H n regions remains a puz- 
zle. We have discussed many effects that may affect the determi- 
nation of the neutral gas abundances such as the source extent, 
ionization corrections, depletion on dust grains, stellar contam- 
ination, and line saturation. The underabundance of a factor 10 
found in the H i region of BCDs is however so dramatic that 
it seems difficult all these effects can make us miss the detec- 
tion of 90% of the metals. Therefore, this offset appears to be 
real and we discuss in the following possible physical expla- 
nations. While self-enrichment of the H n regions, as proposed 
by Kunth & Sargent (1986), naturally explains why metals are 
more abundant in the ionized gas, this scenario is difficult to rec- 
oncile with the homogeneous abundances observed in discon- 
nected H n regions of most dwarf galaxies (e.g., Dufour 1986; 
Russell & Dopita 1990; Skillman & Kennicut 1993; Kobulnicky 
& Skillman 1996; Kobulnicky & Skillman 1997; Legrand et 
al. 2000; Noeske et al. 2000). In fact, the observed dispersion 
is probably dominated by variations of the electron tempera- 
ture and not by intrinsic chemical inhomogeneities. Much less 
frequent is the observation of local discontinuities, with only 
NGC5253 (Welch 1970; Walsh & Roy 1989; Kobulnicky et 
al. 1997; Lopez-Sanchez et al. 2007), IIZw40 (Walsh & Roy 
1993), IC 4662 and ESO 245-G05 (Hidalgo-Gamez et al. 2001). 
Still, the local pollution could be generally explained by Wolf- 
Rayet winds (Brinchman et al. 2008) rather than by SNe explo- 
sions. 

Based on the gas kinematics in the dwarf irregular 
NGC1705, Heckman et al. (2001) suggested that the absorp- 
tion of neutral species could arise in small dense clumps evolv- 
ing in a superbubble interior. The swept-up clumps have likely 
the same metallicity as the ionized gas of the H n region, but 
more metal-poor gas could lie further away, diluting the inte- 
grated abundances. Alternatively, Cannon et al. (2005) proposed 
that the neutral gas probed by FUV absorption lines comes from 
a low abundance H i halo, while metal-rich ionized gas lies close 
to the stellar clusters. The metallicity offset observed between 
the neutral gas and the ionized gas seems to be independent of 
the chemical age of the galaxy, so that abundances have to be 
kept high in the star-forming regions and low in the halo over 
extended periods. This implies that star-formation has to remain 
localized over the galaxy's history. 

Our results and calculations (Sect.|6]l suggest that most of 
the metals released by a starburst episode could mix with the H i 
region, only increasing slightly its metallicity, while a small frac- 
tion (on the order of 1 %) mixes locally within the star-forming 
regions. It must be stressed that these metals were produced by 
previous bursts. In such a scenario, low-level star-formation is 
not required to explain the metallicity of the neutral phase (ex- 
cept if galactic winds allow most of the metals to escape in the 
intergalactic medium). Is this scenario consistent with the ob- 
servation of homogeneous abundances in H n regions of a given 
star-forming galaxy? The fact that star-forming episodes occur 
simultaneously indicates that these episodes were probably trig- 



gered together, or that they triggered each other over several 
lOMyr. Hence it is in fact reasonable to expect relatively sim- 
ilar abundances. 

The contrasting result of Bowen et al. (2005) in the dwarf 
spiral galaxy SBS 1543+593 poses a fundamental question: is 
the metal deficiency of the H i gas specific to BCDs? Although 
Lebouteiller et al. (2006) found that the neutral phase toward 
NGC604 in the spiral galaxy M33 is more metal-poor than the 
ionized gas of the H n region, i.e., being consistent with results 
in BCDs, the situation might be different because of the pres- 
ence of high-velocity clouds along the lines of sight. Rotation 
in spirals provide an efficient mixing throughout the galaxy. 
Furthermore, star-formation keeps occurring along the spiral 
arms. In opposition, BCDs have weak velocity field, and star- 
formation usually occurs in one or several knots (e.g., Papaderos 
et al. 1996). The different mode of star-formation could be re- 
sponsible for substantial differences in the dispersal and mixing 
mechanisms. It seems necessary to study the metal abundance in 
the neutral gas of dwarf irregular galaxies which are not BCDs. 
Studies of other star-forming objects thus ought to be performed, 
which will be possible thanks to the interactive database tool de- 
signed by Desert et al. (2006a; 2006b). Analysis is currently in 
progress. Besides, one of the most complication in the studies 
of absorption-line spectra of galaxies is the presence of unseen 
components, possibly saturated. This caveat is usually solved 
by the use of the weakest observed lines. In that vein, the up- 
coming installation of the Cosmic Origin Spectrograph onboard 
the Hubble Space Telescope will allow the sulfur abundance de- 
termination in the neutral phase through the interestingly weak 

5 ii /U256 multiplet (Pettini & Lipman 1995). The oxygen abun- 
dance will be measured with a better accuracy than with FUSE 
thanks to the O i /U356 line (Meyer et al. 1998). Finally, the ex- 
amination of phosphorus lines will confirm its the possible use 
of phosphorus abundance as a metallicity tracer as suggested by 
Lebouteiller et al. (2005). This is particularly important for the 
analysis of the FUSE database since O i lines are generally sat- 
urated or close to saturation. Finally, it must be stressed that the 
ideal case study for deriving abundances in the neutral gas of 
BCDs is to observe the ISM toward a quasar line of sight (Kunth 

6 Sargent 1986). The method was successfully used in the dwarf 
spiral galaxy SBS 1543+593 by Bowen et al. (2005) but unfortu- 
nately, no alignment BCD-quasar has been observed so far. The 
study of such an alignment would cancel the selection effect in 
FUSE spectra to probe gas toward star-forming regions within 
the galaxy. 

8. Summary and conclusions 

We have analyzed the FUV FUSE absorption spectrum of the 
BCD Pox 36. The following results are found: 

- We detect the interstellar absorption lines of H i, N i, O i, Si n, 
P n, Ar i, and Fe n. Pox 36 appears to be have a particularly 
low foreground gas content, which allows us to measure col- 
umn densities with minimal systematic errors due to hidden 
saturation effects. 

- Metal column densities have been derived from the line pro- 
files. The interstellar H i lines are strongly contaminated by 
stellar absorption. We have modeled the spectral continuum 
using TLUSTY models. The best model is constrained by the 
high-order Lyman line Lye and gives a photospheric temper- 
ature of w 30 000 K (B0 type). After correcting Lyfi for stel- 
lar absorption, we have used the interstellar damping wings 
to infer a H i column density equal to 20.3 ± 0.4. 
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- The observed FUV continuum level is consistent with ~ 300 
BO stars. This represents only the unextincted stars. 

- Abundances have been computed from the column densities 
using the dominant ionization stage. The large [O/Fe] ratio 
in the ionized gas of Pox 36 is interpreted as caused by iron 
depletion onto dust grains. In contrast, the neutral gas does 
not show any depletion, which is attributed to the low dust 
content along the observed lines of sight. 

- We find that the neutral gas is metal-underabundant by a fac- 
tor 5 to 8 as compared to the ionized gas of the H n regions. 
This corresponds to a metallicity of « 1 /35 Z Q in the neutral 
gas as compared to ~ 1/5 Z© in the ionized gas and imply 
that the neutral gas is less chemically processed. 

- It is shown that the metal deficiency in the neutral gas can 
be achieved and maintained if most of the metals released by 
star-formation episodes mix with the H i gas, increasing only 
slightly its metallicity, while a small fraction (1%) could mix 
locally and enrich greatly the ionized gas. The metals should 
come from previous bursts which occurred around the same 
location as the present one. 

- Based on the sample of the seven BCDs studied by FUSE, we 
have confirmed that there is a systematic underabundance by 
a factor ~ 10 of the neutral gas as compared to the ionized 
gas in the H n regions. Besides, [Ar/H] and [N/H] seem to 
reach a minimum value for the lowest-metallicity objects (~ 
1 /8 Z ), which could be due to a default enrichment. 

- There is no convincing evidence that the neutral gas in BCDs 
contains pristine gas pockets, i.e., completely metal-free. 
However, it could be that the neutral gas of BCDs has not 
been affected yet by the star-formation episodes, and that its 
metallicity is that of the IGM. 
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Appendix A: Synthetic H i profiles 

We built H i profiles of Lye, ft, y, and 6 by summing the optical 
depths of the stellar models and of synthetic interstellar lines. 
Profiles were calculated for various column densities and allow 
one to know, depending on the stellar temperature, which line 
can be safely used to constrain the interstellar H i column den- 
sity. From the results of Figs. |A.l| and |A.2| it is clear that: 

- Lya can be used for any stellar temperature, 

- Ly/3 can be used for T > 30 000 K, 

- Lyy and Ly<5 do not provide any constraint, even for high 
stellar temperatures. This is because these lines are saturated 
without any damping wings. 
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Fig. A.l. Prediction of the total H i profile (stellar and interstellar) of the Lyman lines. All stellar models are for 1/1000 Z Q metallicity 
except for the model at 22 000 K for which we take 1/5 Z Q . 
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A.2. See Fig. lA.ll for the plot description. 



